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ABSTRACT 

Aims. We present the analysis of a large sample of gamma-ray burst (GRB) X-ray light curves in the rest frame to characterise their 
intrinsic properties in the context of different theoretical scenarios. 

Methods. We determine the morphology, time scales, and energetics of 64 long GRBs observed by 5w(/f/XRT without flaring activity. 
We furthermore provide a one-to-one comparison to the properties of GRBs witli X-ray flares. 

Results. We find that the steep decay morphology and its connection with X-ray flares favour a scenario in which a central engine 
origin. We show that this scenario can also account for the shallow decay phase, provided that the GRB progenitor star has a self- 
similar structure with a constant envelope-to-core mass ratio ~ 0.02 - 0.03. However, difficulties arise for very long duration (/,, > 10"* 
s) shallow phases. Alternatively, a spinning-down magnetar whose emitted power refreshes the forward shock can quantitatively 
account for the shallow decay properties. In particular we demonstrate that this model can account for the plateau luminosity vs. end 
time anticorrelation. 

Key words, gamma-ray: bursts - radiation mechanism: non-thermal - X-rays 



1. Introduction 



Since the launch of the Swift satellite in 2004 dGehrels et alj 
l200l . the evolution of the gamma-ray burst (GRB) X-ray 
light curves has appeared to be more complex th an previously 
thou ght. The different shapes of the light curves dNousek et alj 
I2OO61) and the existence of X-ray flares superimposed on the 
smoot h continuum (Chincarini et al. 2007, 2010) up to very late 
times dCurran et al.ll2008l : lBernardini et alJI201 1) have prompted 
a large number of studies conducted on different samples of 
X-ray light c urves in order to sin gle out their main properties 
('Zhang et al."2006"; 'O'Brien et al. 2006; 'Willingale et al."2007t 
.Zhang et al. 2007; Liang et al. 2007, 2008, 2009; Gehrels et ID 
r2008l:lRacusin et alJl2009l: lEvans et al.ll2009l) . 

Analysis of the observational properties of 
flares ([C hincarini et al. 2007, 2010; Margutti e t all 
iBernardini et all 120111; |Margutti et al. 2011) allowed us 
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establish a connection with the prompt emission. In particular, 
the extension of the prompt lag-luminosity relation to X-ray 
flares strongly suggests that X-ray flares and prompt emission 
puls es are produced by the same mechanism (Margutti et al. 
I2OIQ) . The nature of the underlying X-ray continuum is still 
debated. It was immediatel y clear soon after th e first set of 
observations by SwiftlXKY dBurrows et al. I l2005h and redshift 
measurements that the light curves differed from a simple 
power-law behaviour an d that some of the characteri stics repeat 
in a systematic way (see iNousek et al . 2006; Zhan g et al.ll2006t 
lO'Brien et al]l2006l) . After six years of data collection, we now 



have a much larger sample for investigating the morphology of 
the X-ray light curves, and a rather large subsample of redshift 
measurements for this analysis in the rest frame. The rest frame 
allows us to understand the intrinsic properties, the timescales 
involved, and the energetics. In fact, the parameters derived 
are directly related to the physics of the phenomenon, and the 
distributions are not affected by the distribution of redshift. 
This is the approach we adopt in this work. We select light 
curves observed by XRT without flares, in order to avoid any 
kind of contamination from other emission components. When 
necessary we c ompar e our results with the sample analysed in 
iMargutti et all d2011') of GRBs with flares, providing for the 
first time a one-to-one comparison of light curves with and 
without flares. 

The initial steep decay of the afterglow is commonly due 
to the tail of the prompt emission (Kumar & Panaitescu 200C|), 
and the plateau phase is often connected with external shock 
emission with an injection of energy. However, the standard 
mo del for GRBs s uffers from several problems (for a review 
see Lvuti kovll2009l) . Here we aim at exploring alternative mech- 
anisms to shape the light curve, such as the accretion of the 
stellar material left behind by the collapse of the massive star 
(iKumar et alJ l2008|)_ or the rotat ional energy of a newly-born 
magnetar (see e.g. lMetzger etalj 2010, and references therein). 
We show indeed that an accretion model would explain the X- 
ray light curves well enough, as well as the properties of flares 
as discussed by IMargutti et al.l d201 Ih . The plateau phase agrees 
with the assumption that new energy is injected from a spinning 
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down neutron star into the forward shock. In particular within 
this framework we demonstrate that this model can account for 
the anticorrel ation between t he plateau l uminosity and its end 
time found bv iPainotti et alJ (l2008ll2()T()h . 

In Sect. |2] we describe the light curve extraction and the lu- 
minosity calibration procedures. In Sect.|3]we describe the sam- 
ple selection criteria, the fitting procedure and the main results 
about the morphology and spectral properties of the sample. In 
Sect. |4]we present the analysis of the energy of the sample com- 
pared with the prompt emission energy. In Sect. |5] we compare 
the present sam ple with the light curves with flares analysed in 
iMargutti et al.l (1201 1). In Sect. |6]we discuss our findings and in 
Sect. |2]we draw our main conclusions. We adopt standard val- 
ues of the cosmological parameters: Ho - 70 km s"' Mpc"', 
Qm = 0.27, and Qa = 0.73. Errors are given at 1 cr confidence 
level unless otherwise stated. Times are in rest frame unless oth- 
erwise stated. 

2. Data reduction 

SwiftfKKl data were processed with the latest version of the 
HEASOFT package available at the time of writing (v. 6.9) and 
corresponding calibration files: standard filtering and screening 
criteria were applied. SwiftlXKY is designed to acquire data us- 
ing different observing modes depending on source count rates 
to minimise the presence of pile-up. However, we applied pile- 
up corrections when nece ssary, following the prescriptions by 
iRomano etal] (12006') and ' Vaughan et"aLr(l2006i) . We extracted 
count-rate light curves in the total XRT 0.3-10 keV energy band. 
Background-subtracted, vignetting corrected light curves were 
then rebinned so as to assure a minimum signal-to-noise equals 
to four When single-orbit data were not able to fulfil the signal- 
to-noise requirement, data coming from different orbits were 
merged to build a unique data point. Count-rate light curves were 
then converted into luminosity curves using the spectral infor- 
mation derived from a time-resolved spectral analysis where the 
spectral evolution of the s ource, if present, is pr operly accounted 
for. We refer the reader to Margu tti et alJ (1201 Ol) for any detail on 
the light curves and spectra extraction. 

3. Temporal and spectral properties of the X-ray 
light curves 

The sample consists of 64 long GRBs observed by SwiftlXKY 
from April 2005 to April 2010. Among all the Swift GRBs ob- 
served in this period, we selected: 

- GRBs with redshift measurements, 

- GRBs without flaring activity in their X-ray afterglow. 

For a list of the GRBs of the sample see Table [T] 

3. 1. Fitting procedure and classification of fhe liglit curves 

The luminosity light curves in the rest frame were fitted adopt- 
ing a simple power-law or a smoothly-joined broken power- 
law model (see Appendix lAlPt a piece of software automati- 
cally identifies the smooth behavio ur of t he X-ray light curve 
of GRBs. We refer to Margutti et aD (1201 ih for details on the fit- 
ting procedure. To reduce the number of free parameters we fix 
the values of the two smoothing parameters di and d2- 

' The reported temporal index of the third segment of the smoothly 
double-broken power law is not the best-fitting value but the asymptotic 
slope (see Appendix lAt. 



Table 1. List of the 64 GRBs of the sample, redshift, and 
morphological type. The GRBs of the golden sample are in 
boldface. 

GRB z' Type GRB Type 



for the redshift measurement, we refer to the values reported in 
the GRB Circular Notices. 

Table 2. Median values and standard deviations (cr) of the 
temporal slopes and break times for the different morpholog- 
ical types (see Fig.[T]i. The values in round brackets refer to 
the golden subsample. 





Type 


Type la 


Type lb 


Type II 




1.29(1.28) 


0.68 (0.80) 


2.9 (-) 


3.3(3.3) 


O-al 


0.27 (0.30) 


0.27 (0.38) 


1.3 (-) 


1.2(0.7) 


(02) 




1.37(1.44) 


0.77 (-) 


0.43 (0.57) 


0-a2 




0.37 (0.27) 


0.37 (-) 


0.32(0.27) 


(as) 








1.5(1.33) 










0.4 (0.23) 






3.25 (3.08) 


2.28 (-) 


1.95(1.85) 






0.73 (0.74) 


0.64 (-) 


0.45 (0.16) 


{Log[tt2/s]) 








3.78 (3.56) 


O" L„s[tb2/s] 








0.50 (0.31) 



The observed light curves can be classified on the basis of 
the best-fitting functior0 (see Fig. [TJ: Type shows a single 
power-law decay; Types la and lb show one break in the light 
curve and a shallow-to-steep or steep-to-shallow transition, re- 



- The full list of parameters is available by requesting it to 
grazia.bernardini @brera.inaf .it 
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Fig. 1. Examples of Type (GRB080804), Type la (GRB081222), Type lb (GRB090814A), and Type II (GRB061021) light curves 
and the respective best-fitting functions and parameters. 



spectively; Type II shows two breaks in the light curv^l A three- 
break model is not required to improve the best fit of any light 
curve of the sample. The present sample contains 22 Type light 
curves, 20 Type la, 7 Type lb, and 15 Type II. In Table|2]we list 
the median value and standard deviation of the temporal indices 
of the power-law segments and the break times of the different 
types. 

For a proper classification of the light curves, we define a 
golden sample composed of 16 GRBs that have been observed 
by XRT from tsian < 40 s to f ^o^ > 10^ s in the rest frame (bold 
in Table [T]|. The limiting start (stop) time has been chosen to be 
much shorter (longer) than the median first (second) break time 
of the Type II light curves (see Table|2]i. 

The presence in the golden sample of five Type la light 
curves suggests that they form a distinct class and are not simply 
Type II observed after the first break. In fact, in those cases the 
observations start time is > 2 cr, so earlier than the first break 
time of Type II, which corresponds to a probability of ~ 5%. 
Therefore, we expect < 1 Ught curve in our Type II sample that 
can be included in the Type la golden sample. A K-S test ap- 
plied to the temporal indices a'" (the first, shallow segment) of 
Type la and a!^ (the plateau) of the Type II for the entire sample 
does not favour a different parent population between the two (P- 
value= 6%). Another possibility is that Type la light curves are 
Type with a jet break. Although we cannot exclude this in some 
cases, the K-S test comparing the temporal indices a'" of Type 
la and of the Type for the entire sample gives a significant 
difference between the two (P-value= 2.2 x 10"*). Moreover, the 

^ T his classification, analogous to the one presented in IE vans et alj 
( l2009ll . corresponds to types d,c,b,a, respectively. 



median break time of Type la ^f^ly - 1795 s is one order of 
magnitude lower than the me dian value fo r jet break times of 
pre- and post-Swift afterglows (iRacusin et al...2009.) . 



The golden sample does not include any Type lb light curve. 
However, in three cases we have observations up to tsiop > 10^ 
s, thus implying that if they are Type II for which we missed 
the second break, this occurs remarkably later than in the Type 
II light curves of the sample (f^j = 6026 s). The K-S test com- 
paring the temporal indices a'2 (the second, shallow segment) 
of Type lb and a!^ (the plateau) of Type II for the entire sample 
gives a 49% probability that they represent the same population. 
The K-S test comparing the temporal indices a['' (the first, steep 
segment) of Type lb and a[' (the steep decay) of Type II for the 
entire sample gives a probability of 99% that they represent the 
same population. A connection is also found between the first 
break time of Type lb t^^ and of Type II tl\ (P-value= 63%). 
Therefore, in what follows we refer to the ai index of both 
Type lb and Type II as the steep decay of the GRB light curves. 
Likewise, we need to further investigate the nature of the shal- 
low decay of Types lb a'^ and la a'" with respect to the plateau 
phase a2 of Type II. 



Also the temporal indices of Type 0, aj" of Type la, and 
a'J of Type II show significant similarities (see Table |3] for a 
summary of the relative K-S test probabilities). We therefore 
group these segments as the normal decay of the GRB light 
curves. 
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Table 3. P-value of the K-S test comparing the temporal indices and break of different types of light curves in the sample. 
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3.2. Spectral properties of the light curves 

The photon indices F have been obtained by extracting high- 
quality spectra within several subsequent time intervals under 
the condition that 2000 photons are included in each time inter- 
val. Then, these spectra were fit by an absorbed power law with 
variable Nh- The behaviour of the spectral index y6 = F - 1 as 
a function of time for each Type of light curve is portrayed in 
Fig.|2] We observe that for Type light curves the spectral index 
is almost unchanged for the full duration of the event apart from 
an initial evolution. The same situation occurs for Type la light 
curves. Type lb shows variability in the spectral indices prior to 
the first break time. This behaviour is analogous to Type II light 
curves, which show variability before the first break time and no 
significant variations before and after the second break times. 

The similar behaviour of the spectral indices of the steep de- 
cays of Types lb and II enforces their association, established in 
the previous section. The other segments do not show any sig- 
nificant spectral evolution, and they are all consistent with ^ ~ 1 
dButler 2007). Type II fight curves tend to converge to a mucfi 
narrower distribution of p than in other cases. In what follows, 
we consider the mean value of the spectral indices of all tfie spec- 
tra covering tfie same time interval as tfie spectral index of eacfi 
segment of eacfi ligfit curve. Tfiis is an excellent approximation 
in all cases but tfie steep decays, wfiere tfie average value is not 
always representative of tfie observed trend. 

4. The energy output 

Tfie energy output of tfie different types of X-ray ligfit curves 
fias been calculated as tfie integral of tfie best-fitting function 
over tfie obse rving time, and isotropic emission is assumed. A K- 
correction ( Bloom et al.ll2001ft fias been applied in order to eval- 
uate tfie energies in a common rest frame energy band (1.9-11.2 
keV). For th e prompt er nission isotropic ene rgy Eiso, we refer to 
lAmati et all (12008.) and lKrimm et all (l2009l) . and to lNavaet all 



(1201 oh for those GRBs observed by Fermi/GBM whose prompt 
emission spectrum has been fit with a Band function. Figure[3]a 
shows the distribution of the total energy output of the X-ray 
light curves Ej^ as a function of Etso- It shows a moderate, pos- 
itive correlation (with a Spearman rank p = 0.56, see also 
IWiUingale et"ani200l . 

We overplotted in Fig. [3] a the total fluence in the XRT energy 
band (0.3 - 10 keV) as a function of the prompt emission fluence 
detected by BAT (15-150 keV) for the complete sample of long 
GRBs observed by XRT until December 2010 without flares and 
with a complete monitoring of the X-ray light curve from before 
300 s in the observer frame until it reaches the background level. 
These data are provided by the XRT Catalog ue (Margutti et al., 
in preparation; see also IZaninoni et all [201 lb . The existence of 
a positive correlation between the X-ray and prompt emission is 
also apparent in the observer frame (p = 0.51). 

The X-ray energy that we computed for the present sample 
is a lower limit of the total, since we only consider the emission 
during the observing time in a limited energy band. To partially 
overcome this limit, we consider the light curves belonging to 
the golden sample separately, and we integrate those curves be- 
tween the end of the prompt emission (the rest-frame Tgo) and 
10^ s. The positive correlation between the prompt emission and 
the X-ray energy for this subsample is enhanced (p - 0.59, see 
Fig.Ela). 

We now consider each light curve type separately. Type 
light curve energy Epi shows a moderate correlation with the 
prompt emission energy, but the paucity of Type ligfit curves in 
tfie golden sample does not allow us to be more conclusive (see 
Fig.[3]b). Tfie same situation is found for tfie steep decay of Type 
II and lb energy Esieep (see Fig.|3]d) and for tfie normal decay of 
Type II and la energy E„orm (see Fig.|3]f). Wfien we consider tfie 
corresponding fluences as a function of tfie 15-150 keV fluence 
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of the prompt emission for the sample defined abovfl we do not 
find any trend toward Type fight curves (p = 0.35) or toward 
the normal decay (p - 0.33), and a moderate correlation for the 
steep decay (p - 0.49). 

The energy output of the plateau phase Epia,eau in Type II 
light curves is positively correlated with (see Fig. [3]e). In 
this case we can consider the entire sample (the initial and end 
time of the plateau are always known), and we find a linear 
dependence: Epia,em, with Epia,eau/Ei,o ~ 2%, which 

is roughly the s ame order of magn itude than E gnrel Eis,-, fo r 
early-time flares jChincarini et alj|201 0; Berna rdini et al 
A correlation is also found between the prompt emission en- 
ergy and the initial shallow decay of Type la energy E shallow 
(see Fig. [3]c). If we consider only the golden sample, we find 
that this dependence is almost linear: Eshaiiow °^ Eiso, with 
Eghaiiow/Eiso ~ 3%. Both these correlations are also apparent in 
the corresponding fluence-fluence plot (p = 0.58 and p - 0.56, 
respectively). We have to caution that a prompt monitoring of 
the light curve is necessary in order to distinguish Types la and 
II. Therefore, the fluence-fluence plot for the shallow decay can 
contain spurious events due to a misclassification. 



The light curves of the XRT Catalogue are fitted and classified with 
the same procedure as adopted in this work, therefore w e are comparing 
equiv alent quantities for each part of the light curves (see lZaninoni et all 
l20Tlh . 



5. X-ray light curves with and without flares 

The selection of X-ray light curves without flares allows us to 
investigate whether the X-ray continuum knows about the pres- 
ence of flares. For this, we compare our result s with the cor- 
responding values of the sample considered in iJVIargutti et al.l 
dlOl ih of 44 GRBs that exhibit flaring activity. The major ad- 
vantage is that the fitting procedure and the functions used for 
the continuum underlying the X-ray light curve of that sample 
are the sam^B- 

If we classify the light curves with flares (FLC) on the basis 
of the ir best-fitting function as in Sect. l3.1l we find that the sam- 
ple of iMarputti et al.l ( 1201 l\i contains a majority of Type II light 
curves (72%), 14% of Type 0, and 7% of Types la and lb. This 
repartition is different from what we found in our sample (34% 
of Type 0, 31% of Type la, 1 1 % of Type lb, and only 23% of 
Type II, see also lMargutti et al.ll20ldh . 

We produced a population of "fake" flares on the basis of 
the relations found in Chincarini et al. (2010) (w - 0.2 tpk + 10; 
tdec ~ 2w; Lpk oc r^.v. E^i^^^^ ^ iqSI gj.g) ^jjj^ 4Q ^ 

< tpii < 500 s, and we superimposed these flares to the light 
curves of our golden sample. We find that a "standard" flar- 
ing activity would have been detected above 1 cr in 80% of the 
golden sample light curves. In particular, these flares could have 
been detected in all the five Type la light curves for f > 100 
s. Moreover, if we compute the average luminosity of the light 



' We recall, however, that for the flare sample the luminosity is in the 
common rest-frame energy band 2.2 - 14.4 keV. 
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Fig. 3. a) Total energy of the X-ray light curve vs. prompt emission energy Eiso- The black line marks the condition E^ = Ejso- 
b) total energy of the Type light curves Epi vs. Ej^o- c) total energy of the first, shallow segment of Type la light curves EsiwiUm- 
vs. Eiso- The red line is the best-fit of the golden sample data: Log[E shallow] - (20 + 85) -i- (0.6 + 1 .6)Log[Eiso\- d) total energy of the 
steep decay of Type lb and II light curves Esteep vs. Eiso- e) total energy of the plateau of Type II light curves E plateau vs. Eiso- The 
blue line is the best fit of the complete sample data: Log[Epia,eau] = (5.8 + 1.4) -H (0.85 + 0.03)Log[Eis„]- f) total energy of the final 
normal decay of Type la and Type II light curves E„orm vs. Eiso- The grey dots refer to the fluence-fluence plot of the 278 long GRB 
light curves with (stars) and without (squares) flares from the XRT Catalogue by Margutti et al. (in preparation). 



curves of the golden sample, we find that it has a shallower de- 
cay {{L} oc r^) but is even dimmer than the average continuum 
obtained in Margutti et al.l (1201 li) for t < 200 s. This excludes 
a luminosity bias preventing detection of early time flares in the 
light curves of our golden sample. 

In Fig. |4] we compare the distribution of the steep decay 
power-law index of Type II and lb light curves of the present 
sample (LC) with the corresponding value found for the FLC. 
We did not include those cases in which the prompt emission 



still observed in the X-rays (as GRB060510B). The FLC power- 
law index distribution has a median {a'^^'^ - 2.7 {cr^ - 1.3), 

while the LC sample has {a''^^) = 3.2 (a-^ = 1 .2). The K-S test 
comparing the temporal indices o-^'- and gi ves a probabil- 
ity of 9 % that they belong to the same population. iMar^utti et al.l 
(1201 Ih show that light curves with multiple flares tend to be shal- 
lower than those with a single flare. We then select the FLC with 
a single flare that are expected to be more similar- to LCs and 
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Fig. 4. Probability distribution of the steep decay power-law in- 
dex asleep for the LC (blue lines) compared with the FLC sample 
(red Une). 



compare their power-law indices with the LC distribution. The 
two distributions are more similar, with a median of the single 
FLC distribution (a'l'f) = 3.0 (o-^.f;f = 1.0). The K-S test com- 
paring the temporal indices a^^ and ff^f^^ results in an increasing 
probability that they belong to the same population (27%). 

We aim at investigating now the late-time behaviour of the 
entire sample. For this purpose we computed the average light 
curve for the entke LC sample. After 1000 s we observe that 
(L) oc f a trend similar to the detection threshold fo und 
for the average flare light curve in 'Margutti et al.l (1201 ih . If 
we compare the average luminosity of the LC sample with the 
one of the continuum u nderlying the late-time flares (LFLC) of 
iBernardini et al.l (1201 ll) . we find that the temporal behaviour is 
very similar but the present sample is brighter, with L^^^'^ j L^c" ^ 
20%. 

We turn now to the energy output of the X-ray light curves 
with flares. We restrict our analysis to a flare golden sample 
(FGS) of 17 light curves, defined with the same criterion as was 
introduced in Sect.|3] and we calculate the energy of the underly- 
ing continuum in the common rest frame energy band (1.9-11.2 
keV) as in Sect.|4] We find that the distributions of both the X-ray 
energy and the isotropic prompt emission energy are similar for 
the two samples. A correlation between the and Eiso is present 
{p - 0.70, see Fig.O, as well as between the total fluence in the 
XRT energy band (0.3 - 10 keV) and the prompt emission flu- 
ence detected by BAT (15 - 150 keV) for the complete sample 
of long GRBs with and without flares from the XRT Catalogue 
by Margutti et al. (in preparation) with a complete light curve 
in the sense described in the previous section (p = 0.51). The 
correlation between E plateau and Ej^o is enhanced including the 
FGS data (see Fig.[3]e). The results of previous section are con- 
firmed by adding the light curves with flares also in the fluence- 
fluence plot for the plateau phase (p = 0.50) and the shallow 
decay (p - 0.63). The other fluence-fluence relations portrayed 
in Fig.[3]are not improved significantly with the exception of the 
Type versus prompt emission (p = 0.50). 

6. Discussion 

The analysis of the rest-frame X-ray light curves of our sample 
shows that they can be grouped into four morphological types. 
We found that the steep decay of Types lb and II can be grouped 



together, as well, as the shallow decays in Types la, lb, and II, 
and normal decays in Types 0, la, and II. Spectral evolution is 
associated with the steep decay, while no spectral evolution is 
found during shallow and normal decay. We demonstrated that 

- a quasi-linear correlation exists between the total isotropic 
energy emitted during the prompt phase Ejso and the energy 
released during the X-ray shallow decay phase of GRBs with 
and without flares. This correlation is also apparent in the 
observer frame of a larger sample; 

- the predominance of Type II with flares with respect to the 
present sample indicates that the probability of flare occur- 
rence is different depending on the light curve morphology; 

- the energetic of the X-ray continuum with and without flares 
is similar. Light curves without flares tend to be steeper than 
those with flares. 

6. 1 . The steep decay 

The widespread consensus is that the initial steep decay ob- 
served in Types lb and II is the tail of the prompt emission. 
If the central engine stops its activity at the end of the prompt 
emission, the decline in the X-ray l ight curve is imposed by th e 
high-latitude emission (HLE, r^-^: lKumar & Panaitescull2000h . 
This rapid decline is possible only if Oj > F \ where Oj is the 
jet opening angle and F the Lorentz gamma factor of the ejecta. 
If Oj < F"', the light curve morphology is determined by the 
emission of ejecta that cool vi a adiabatic expansion (AE) (AE; 
iBarniol Duran & Kumaijl20()9|). Variations in the decay rate are 
expected for structured iets dzhang & Meszarosll2002l : iLvutikovl 

2006) . 

In Fig.|5]we show that 32% of the sample agrees with HLE 
within Icr, 32% with AE, and 79% are within the ranges im- 
posed by file structured jet. GRB090618, GRB081203A, and 
GRB081118 are not compatible with any case. If we include 
Type l b and II GRBs of the sample analysed in Margutti et aT| 
( 1201 Ih . we find that 82% are consistent with one of the above 
cases. 

The simplest scenario predicts that no spectral evolu- 
tion should be observed during this stage(Zhang et al. 2007t 
Barniol Duran & Kumar 2009; see, however, Butler & Kocevskj 

2007) . Figure|2]clearly shows that this requirement is never met 
for the Type II and lb light curves of our sample. GRB090618 
and GRB081203A are clear examples: we show in Fig.|5]the dif- 
ferent positions that GRBs assume during the steep decay in the 
o:steep - P Sleep plane due to the spectral evolution. 

The spectral e volution of the early steep deca y and other ar- 
guments (see e.g. ^ Barniol Duran & Kumaill2009 ^ lead us to be- 
lieve that the observed early X-ray decay of GRBs is produced 
by the rapidly declining continued activity of the central engine, 
if we assume that there is a one-to-one correspondence between 
the temporal beh aviour of the central engine act ivity and the ob- 
served emission ('B arniol Duran & Kumad2009h . This is the case 
for a continuous accretion onto the central object, as proposed 
bv lKumar et alJ (l2008l) . The decay indices {{asteep) - 3.15) and 
spectral evolution we observe in our sample support this hypoth- 
esis. In fact, the predicted light curves decay as L oc or steeper 
when the density in the outermost layer of the progenitor star has 
a sharp density profile, or L ~ for a hi gher rotation rate of the 
pro genitor star core (see Figs. 4 and 5 in lKumar et al.ll2008l see 
also lLindner et al ] l2010l for the results of hydrodynamic simula- 
tions). 

Within this framework, flares could be powered by insta- 
bilities affecting the physical source of energy, which creates 
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Fig. 5. Temporal decay index Osieep vs. spectral index /Ssteep of 
the first segments of Type II (green stars) and Type lb (pink tri- 
angles) light curves of the present sample. Grey symbols refer to 
the sample of lMargutti et al.l (1201 1). The black stars track the po- 
sitions for the different y8(f) of GRB090618 and GRB081203A. 
The grey solid line corresponds to the expectation for the HLE: 
a = 2+/3. The shaded area corresponds to the expectation for the 
structured jet (-3 - j0<q'< -1 -yS). The blue lines correspond 
to the expectations for the synchrotron emission from AE ejecta: 
a = 1 .5/3 +1.5 for thick ejecta (solid line) and a = 3/3 + 3 for thin 
ejecta (dashed line). The red lines correspond to the expectations 
for the SSC emission from AE ejecta: a - 1.5/3 + 3.5 for thick 
ejecta (solid line) and a = 3/3 + 5 for thin ejecta (dashed Une). 

the steep decay |Kumar et al.l l2008l see also the results from 
iMargutti et akl |20T H). This is a plausible explanation for the 
higher occurrence probability associated with the steep decay. 
Therefore this scenario also has the merit of accounting for the 
connection steep decay/flaring activity. 

6.2. The plateau phase 

The intermediate shallow decaying phase observed in Type 
II light curves is interprete d as an injection of energy into 
the forward shock (see e.g. IZhang et"an l2006t and references 
therein). The absence of significant spectral evolution during this 
stage agrees with the expectations from forward shock emis- 
sio n (see Fig. IS. Alte rnative possibilities have been discussed 
by iFan & PiranI (|2006|) . including a small fraction of electrons 
being accelerated, evolving shock parameters, low non-constant 
dens ity of the interstellar me dium, and a magnetised outflow (see 
also iLyutikov & Blandfordll2003 ). These authors conclude that 
the most viable scenario is the energy injection model, although 
it cannot self-consistently account for both the X-ray and optical 
afterglows in all cases (see e.g. iPanaitescu & Vestrand.2011,) . 

6.2.1. Energy injection from spinning-down neutron star 

A possible source of energy injection is th e power emitted 
by a spinning-dow n newl y-bor n magnetar (Dai & Lu| Il998t 
IZhang & MeszarosI 120011; ICors i & Meszaros 2009) that re- 
freshes the forward shock (see iDalFOsso et al.l |201 ll for an 
analytic treat ment of this problem). We applied the solution 
of [PairOsso et al.l (1201 Ih to our sample (details are given in 
Appendix iBt. and we find that it fits the observed light curves 



well enough. This is independent of the presence of flares (see 
Table IrTT i. provided that the light curve's normal decay agrees 
with the forward shock model prediction^ 

The major advantage of this interpretation is that all the 
plateau properties are directly related to the central engine and, 
consequently, to the prompt emission. This explains the connec- 
tion between the plateau energy and the prompt emission energy. 
The basic information can be derived from two main quantities: 
the magnetic field B and the period P of the pulsar. In partic- 
ular, they are all that is needed to explain the anticorrelation 
between the plate au luminosity Ln and the second break time 
tb2 - tp found by Dainotti et all (l2008i 1201 0). Starting from the 
distributions of B and P from our sample (see Appendix iBl). we 
obtained the normalisation, the slope, and the scatter of the ob- 
served anticorrelation (Fig.|6l for details see Appendix iBl). We 
underline that these distributions are within the range of val- 
ues expected for newly born, millisecond spinning magnetars 
(.Duncan & Thompson 1992) . 

We plot in Fig. |6] the luminosity at the break time L{thi) of 
Type la light curves in the rest frame energy band corresponding 
to I = 2.29, and we find that they follow the same luminosity- 
time anticorrelation and that they can be interpreted as the spin- 
down luminosity of a milUsecond pulsar for the same distribu- 
tion of B and P as for Type II. We add to Fig. |6]last Type lb 
observation as a lower limit on the end of the injection phase 
and an upper limit on the luminosity. The possibility of having 
injection times up to 10^ s, as observed in Type lb light curves, 
is allowed within reasonable values of the magnetic field and pe- 
riod (see Fig.|6]). However, the upper limit on luminosity found 
for some Type lb is much lower than the expected one, unless we 
assume that the injection time is > 10'' s. We argue that a differ- 
ent beaming factor and/or efficiency in converting the spin-down 
power in X-rays may account for such Type lb light curves. 

The main constraints of the magnetar model is related to the 
energy budget. The maximum energy emitted in such a model is 
a few 10^^ erg and limited by the maximum r otation energy at- 
tainable by a rotating neutron star (IUsovll992h . The energy bud- 
get strongly depends on the uncertain estimate of the jet angle, 
however in a few cases the released energy may be high enough 
to challenge the model dCenko et al.ll2010l) . 

6.2.2. Tine plateau pliase in accretion models 

The plateau phase in the accretion model of iKumar et al.l (l2008h 
is triggered if the progenitor star has a core-envelope structure. 
An analysis of Type II light curves with this model has been 
presented by Cui et al. (2010). They inferred properties of the 
progenitor star based on the timescales and energetics, although 
they assume that the initial steep decay is not produced by the 
central engine. This approach is very effective in obtaining some 
indications about the progenitor star structure. 

In Sect. ID we show that a quasi-linear correlation exists be- 
tween shallow decay energy and the prompt emission total en- 
ergy for the light curves with and without flares, extending the 
relation between the plateau energy and the prompt emission en- 
ergy observed in 15 - 150 keV p reviously found bv iLian^ et al.l 
(.2007) (see also ICui et al.ll2010l) . It implies that the progenitor 
star has a self-similar structure, with a constant envelope-to-core 
mass ratio ~ 0.0 2 + 0.03. This is valid for Type II, as shown by 
Cui et al. (1201 Oh . but also for Type la, where similar injection 



* The applicability of the magnetar model to GRBs displaying a very 
steep decay after the plateau inconsistent with th e forward shock emis- 
sion has been investigated bv lLvons et alJ ( 120101) . 
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Fig. 6. Luminosity at the end time versus end time of the plateau phase. The black squares are the sample analysed bv^D ainotti et alJ 
(l20Toh . The red squares are Type la light curves of our sample and the pink squares are the luminosity of the last observation in Type 
lb light curves in our sample. The grey dots are the 100000 simulations of the luminosity at the spindown time and the spindown 
time assuming that the magnetic field and the NS period are normally distributed around the mean values we found for our sample. 
The blue line marks t he region that includes 99% of the simulations. The green stars are the values found for our best fit with 
iDall'Osso etal] (1201 Ih model. 



times (see Table |2]i point to similar dimensions of the progenitor 
star. The slightly different distribution in the power-law indices 
of the two cases (see Table|2]i can be ascribed to a p ossibly difi'er- 
ent t ypical accretion rate of the envelope (see e.g. iKumar et al.l 
I2OO8 . their Fig. 7). 

However, in Sect. 13.11 we show that the shallow decay in 
Type lb is closely related to the one in Types la and II. This 
poses severe problems to this scenario. In fact, it is very diflicult 
to obtain a shallow decay duration > 10^ s, since such an ex- 
tended envelo pe is more probably ejected during the main burst 
(iKumar et al.li2008.) . 

6.3. The normal decay 

The forward shock emission is expected to give a significant con- 
tribution to the X-ray light curve after the fading of the prompt 
emission (see e.g. Zhang et al. 2006, and references therein). 
The shape of the ligh t curve depends on the density profi le of 
the ambient medium dSari et alJll998HChevalier & LillioOOl) and 
the tempo ral and spectral indices sh ould satisfy specific closure 
relations ("Zhang & Meszarosl l2004l) . Figure [T] shows that 80% 
of Type II light curves matches the expectations of the standard 
afterglow model in the slow cooling regime within Icr, 90% of 
the second segment of Type la, and 77% of Type light curves. 



This connection is favoured in Type la and II light curves by the 
absence of spectral variation during the transition from the pre- 
vious injection phase (see Fig.|2]l.Kwe include Type II GRBs of 
the sample analysed in Margutti et al.l (1201 ll) . we find that 84% 
are consistent with one of the above cases. 

If the origin of the plateau phase is inj ection of energy from a 
millisecond pulsar into the forward shock dPalFOsso et al.l201 iL 
and references therein), after the plateau we expect that the light 
curves follow the standard afterglow scenario, as we find for the 
majority of cases. Alternatively, the accretion model predicts an 
asymptotic behaviour ~ 3^2.7) ^^j-gj. jj^g gj^^ jjjg ^^jj ^^^.j^ of 

the envelope, which covers the range of variabiUty of our a„onnai- 
However, this does not constrain the spectral index of the final 
decay to be related to the one of the previous accretion phase, 
as we found in Fig.|2] On the contrary, since the late time decay 
should be produced by the same mechanism than the early steep 
decay, a similar spectral behaviour should be expected. 

6.4. The occurrence of flares 

The comparison of the properties of the present sample with 
the X-ray light curv es of GRBs showing flaring activity in 
Ma rgutti et al I(|201 ll) shows that flares occur more likely in Type 
II light curves (73% of the flare light curves are Type II, while 
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Fig. 7. Temporal decay index a„onn vs. spectral index /3„orm of 
Type (red squares), of the second segment of Type la (blue 
triangles), and the third segment of Type II (green stars) light 
c urves of the present sa mple. Grey symbols refer to the sample 
of'Marguttie talJ(l201lh . The shaded areas correspond to the ex- 
pectations of the standard afterglow models in the slow cooling 
regime: the lower area only is for ISM model, the intermediate 
area plus the lower area are for the wind model (see e.g. Table 
2 in Zhang et al. 2006.) . the upper area is for the post jet-break 
(ISariet all 19991) 



only 23% of the lig ht curves without flares, as first noticed in 
lMarguttietal]|2010 V We argue that the presence of the steep de- 
cay is related to suitable conditions for the occurrence of early- 
time flares. A connect ion between the fl ares and the steep decay 
has been found in Mar gutti et alj ( 1201 1): the average luminosity 
of flares decays in time as the average slope of the steep decay. 
In that work we also noted that light curves with multiple flares 
are associated with shallower decays. A hint of a steeper distri- 
bution of the power-law index for the light curves without flares 
is found in Fig.H] as if there is a transition from steep to shallow 
decays depending on the number of flares. However, the moder- 
ate difference is not conclusive. 

The association of the steep decay with the flaring activity 
is an argument for a possible contribution of the central engine 
to the steep decay emission. In this case, the slope of the steep 
decay in the accretion model depends on the rotation rate of the 
progenitor cor e: the faster the rot ation, the shallower the decay 
(see Fig. 5 in iKumar et al.l 12008 "). A fast rotation of the infall 
material is needed to create the suitable conditions for insta- 
bihties, which are supposedly respon sible for accretion models 
dPerna et al ] l2006HKumar et al.ll2008l) . A change in the rotation 
rate of the core produces a change in the peak jet power and of 
th e duration o ver which the luminosity is high, as shown in Fig. 5 
of iKumar et al. (2008). However, we do not observe any signifi- 
cant difference in the L^jt and Ejso distributions of the GRBs with 
and without flares. 

In lBemardini et al] (1201 Ih we showed that the decoupling of 
the evolution of the peak luminosity of late-time flares from the 
underlying continuum allows the detection of only the brightest 
flares. Comparing the average luminosity of the present sample 
with the one calculated for the continuum underlying the late- 
time flares analysed in Bernardini et al. (20 11), we observe that 
their temporal behaviour is similar, an d it traces the flare de tec- 
tion threshold luminosity calculated in lMargutti et al.l (1201 ih . 



7. Summary and conclusions 

We analysed 64 long ORB X-ray light curves observed by XRT 
with redshift measurement that do not exhibit flaring activity. 
This allowed us to characterise the morphology and energetics 
of the sample in order to constrain the mechanism that produces 
the X-ray continuum. The light curves of the sample can be di- 
vided into four morphological Types, with similarities between 
different parts. When nece ssary, we compared our sample with 
the one considered in iMargutti et all ( 1201 li) of 44 GRBs that ex- 
hibit flaring activity. We found that 

- a large fraction of the light curves are on average consistent 
with the high latitude emission from a structured jet, but all 
of them show a strong or moderate spectral evolution. The 
accretion model accounts for the observed range of temporal 
decays and the spectral evolution of the steep decay phases 
of Types lb and II; 

- the injection phase observed in Type la and II can be in- 
terpreted as power emitted from spinning-down ultramagne- 
tised neutron star that refreshes the forward shock. We ap- 
plied such a model to our sample and found that it fits the 
data well enough; 

- this scenario accou nts for the o bserved Lp - tj, anticorrela- 
tion Dainotti et al.l (I2008il2010l) : we reproduced the normal- 
isation, slope, and scatter of the anticorrelation starting from 
plausible values of the parameters. The consistency of Type 
la shallow decay with the anticorrelation within 99% con- 
fidence level suggests that they can be reproduced by this 
model with the same distribution of B and P as for Type II; 

- a second phase of accretion of a progenitor with a struc- 
ture core envelope can originate the shallow phase in Types 
la and II. The correlation between the prompt emission to- 
tal isotropic energy with the shallow decay X-ray energy 
implies a self-similar structure for the progenitor star, with 
~ 0.02 - 0.03. However in this scenario, difficulties arise for 
the very long duration {t > 10^ s) of Type lb shallow decay; 

- 70% of the light curves with flares are Type II, while only 
23% of the non-flare sample are, confirming our previous 
finding dMargutti et al.l l2010h . This suggests that Type II 
light curves are related to the suitable conditions for the 
occurrence of flares. The accretion model has the merit of 
proposing a unique mechanism to explain the steep decay 
and the flaring activity; 

The magnetar model is able to explain the different mor- 
phologies of the X-ray continuum and the properties of the 
shallow decay. The main limit in the accretion model is the 
difficulty of achieving very long timescales as observed in 
Type lb shallow decay. The behaviour of the steep decay and 
the suggested connectio n wit h the flaring a ctivit y (see also 
iBarniol Duran & Kumadl2009t iMargutti et al.ll2010l 12011 indi- 
cates that the steep decay may not be viewed as simply the tail 
of the prompt emission. 
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Appendix A: The best-fitting models 

Type light curves correspond to a simple power-law model: 

Lit) = Nr"'. (A.l) 

Type I light curves correspond to a smoothly-joined broken 
power-law model: 



L(t) = N 



t Y' It 
thii \tbi 



(A.2) 



Type II light curves correspond to a smoothly-joined 
double-broken power-law model: 



L(f) = N 











i tf 
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V 



1 + 
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(A.3) 



The a3 reported throughout the paper is not the best-fitting pa- 
rameter ffj that appears in Eq. IA.3I but its asymptotic value 
L{t » thi) «: r'". 



Appendix B: Energy injection from a spinning down 
neutron star 

The injection of energy into the forward shock from a millisec- 
on d spinning, ultramagne tised neutron star (NS) can be written 
as (iDairOssoet al.ll201 ih 



dE ,E(T) IKcJ 

— = L,d{T) - r — — ■ = '— 

dT T (1 + aT)^ 



(B.l) 



where / is the NS moment of inertia, K - B^R^ /{6Ic^) with 
B the dipole magnetic field at the NS pole, R the NS radius, c 
the speed of light, a = 1 /f2 = IKojf is the inverse of the spin- 
down timescale, k' is the radiative efficiency, and E the initial 
blastwave energy. A solution of this equation is (Dal FOsso et alj 
1201 ih 

E{T) = — rr- ^ + — (B.2) 

where Ta is any time chosen as initial condition. The solution 
of the above integral can be expressed in terms of the real val- 
ued hypergeometric function 2Fi(a,b,c;(l + aTY^). The total 
bolometric luminosity of the blastwave is, then 



L{T) = k'E{T)/T. 



(B.3) 



We select 12 Type II light curves from our sample that show 
a well-sampled, prominent plateau (see Table lB.ll ). and we com- 
pute the 0.3 - 100 keV energy band luminosity. We fit this sub 
sample with Eq. lB.3l We use as free parameters B, the NS period 
P = Eo, and k'. We adopt as starting time of the plateau 

phase To the (rest-frame) time in which the spectral index /3 be- 
comes approximately constant. Figure IB. II shows the result of 
the fit in one case, while the best-fit parameters for the entire 
sample can be found in Table IB. II We find that both the mag- 
netic field and the NS period {B = 3.4 x 10'^ G; P = 4.2 ms) 
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Fig. B.l. GRB090618 light curve with the model in Eq.|B3](red 
line) for the best-fit values reported in Table IB. 11 

Table B.l. Table of the Type II GRBs fitted with Eq. iBJl 
z is the redshift, To the initial time, B the best-fit values for 
magnetic field, P for NS period, k' for radiative efficiency and 
Eo for blastwave energy. In the lower panel the same values 
are listed for a set of GRBs with fl ares (see also the f it of 
GRB060729 and GRB061121 in Da H'Osso et all (|20TT1) ). 



GRB 


z 


n(s) 


B(10" G) 


P (ms) 


k' 


£. (10™ erg) 


051016B 


0.9364 


600. 


7.0 ± 1.3 


11.0± 1.0 


0.80 ± 0.36 


0.07 ± 1.3 


05 11 09 A 


2.346 


100. 


2.8 ±0.1 


1.5 ±0.1 


0.42 ± 0.06 


10.0 ± 1.3 


051221 


0.5465 


100. 


4.0 ± 0.0 


15.0 ±0.9 


0.52 ±0.12 


0.3 ±0.2 


060502A 


1.51 


400. 


2.2 ±0.5 


2.7 ± 0.4 


0.25 ±0.13 


3.0 ±0.8 


060605 


3.78 


300. 


4.0 ± 1.8 


2.14 ±0.2 


0.99 ± 0.56 


3.0 ± 1.5 


070306 


1.4959 


1000. 


1.7 ±8.9 


2.8 ± 0.2 


0.99 ±0.33 


1.5 ±0.5 


080707 


1.23 


100. 


5.0 ± 0.8 


10.2 ±0.1 


0.99 ± 0.60 


0.15 ±0.36 


090529 


2.625 


1000. 


1.2 ±0.6 


2.8 ±0.1 


0.89 ±0.17 


2.0 ± 0.8 


090618 


0.54 


160. 


3.0 ± 0.0 


2.5 ± 0.0 


0.81 ±0.01 


12.0 ± 5.0 


090927 


1.37 


1000. 


3.2 ± 0.3 


5.1 ±0.1 


0.99 ±0.12 


0.8 ± 1.4 


100425A 


1.755 


700. 


4.0 ± 0.7 


5.2 ± 0.3 


0.51 ±0.14 


0.9 ± 0.8 


100621 


0.542 


250. 


3.0 ± 0.0 


3.8 ± 0.0 


0.70 ± 0.02 


4.0 ± 0.3 


060526 


3.21 


100. 


3.2 ± 1.3 


5.5 ± 1.0 


0.70 ± 0.06 


1.0 ±0.6 


060714 


2.71 


60. 


5.0 ±2.0 


2.8 ± 0.2 


0.43 ± 0.05 


10.0 ± 6.0 


080310 


2.43 


400. 


6.3 ±0.8 


5.8 ± 1.0 


0.80 ± 0.06 


4.0 ± 1.2 


091029 


3.79 


100. 


1.99 ±0.6 


1.6 ±0.3 


0.32 ±0.01 


1.6 ±0.4 



are within the range of va lues expect ed for newly born, millisec- 
ond spinning magnetars (iDuncan & Thompson 1992] ), and con- 
sistent with the results obtained by .DalFOsso et alj (1201 Ih and 
iLvons etalJ (12010). 

Assuming that the B and P are normally distributed around 
the mean values, we found for our sample B - 3.4 x 10'^ G and 
Log[P/ms] = 0.62, with standard deviations cr^ - 1.5 x 10'^ 
G and criog^pf,,,,] - 0.31. We calculated the spindown timescale 
f2 oc B^/P^ and the corresponding luminosity L(t2) according to 
Eq. lB.3l (rn = 10 s, A; = ^ and Eo - Eo). Using an average spec- 
tral index ^ - 1, we normalised the luminosity in the rest frame 
0.3 - 10 keV energy band c orresponding to the average redshift 
of the iDainotti et alJ (1201 Oh sample, z - 2.29. The results are 
portrayed in Fig. |6l 
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